Soft Spheres Make More Mesophases 
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We use both mean-field methods and numerical simulation to study the phase diagram of classical 
particles interacting with a hard-core and repulsive, soft shoulder. Despite the purely repulsive 
interaction, this system displays a remarkable array of aggregate phases arising from the competition 
between the hard-core and shoulder length scales. In the limit of large shoulder width to core size, 
we argue that this phase diagram has a number of universal features, and classify the set of repulsive 
shoulders that lead to aggregation at high density. Surprisingly, the phase sequence and aggregate 
size adjusts so as to keep almost constant inter-aggregate separation. 

PACS numbers: 61.30.Dk, 61.30.St,61.46.Bc 



Entropy is a potent force in the theory of self- 
assembly. It can be argued, through entropic considera- 
tions alone, that hard spheres will self-assemble into the 
face-centered-cubic (foe) lattice or any of its many vari- 
ants related through stacking faults. As a result, when 
a material exhibits an fee phase, it is often attributed 
to the optimality of the close-packed lattice. When less 
common or less dense lattices are formed, all sundry of 
explanations are invoked, ranging from quantum mechan- 
ics Q , lattice effects , partially filled Landau levels , 
and even soft interactions 0, U While there has been 
concerted effort to tailor the pair interaction to achieve a 
desired periodic arrangement [|| this must be done in the 
context of those packing motifs that arise from generic 
interactions. For instance, it would be no trick to tailor 
a potential to make an fee lattice. 

With this in mind, here we consider a seemingly simple 
extension of the hard sphere model, namely a hard core 
interaction of radius a with a soft shoulder of radius a s 
and height e (HCSS): 
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When <j s /a > 1 this potential models hard spheres with 
a soft pair repulsion and was used to study isostructural 
transitions in Cs and Ce Q. In generic repulsive poten- 
tials, it has been shown that as the range of the soft repul- 
sion grows (corresponding to <r s /a ~ 2) a rich variety of 
density- modulated ground states appear 0,0,01 which 
can be characterized as periodic arrangements of regu- 
lar sized clusters of the original spheres. In this Letter 
we establish a sufficient condition on the pair potential 
for clustering and the subsequent ordering of the clus- 
ters which generalizes results on soft potentials without 



for soft repulsion and use it to study the formation of 
striped phases. We corroborate our analytic treatment 
with numerical solutions that also predict the existence 
of hexagonal and inverted hexagonal phases with both 
fluid and crystalline order in the clusters, as shown in 
Fig. ^ We also present results from Monte Carlo simu- 
lations of the HCSS potential which both stimulate and 
support the more general results. In all cases, we find 
that over the range of stable aggregate structures the 
lattice constant remains fixed while the clusters change 
their size and morphology so as to maintain the average 
intra-cluster density p. 

Why should repulsive potentials lead to aggregation? 
We can understand this in the context of the HCSS 
model: consider a uniform density of spheres at densities 
just large enough so that the soft-shoulders only begin 
to touch the hard cores. The energy can be lowered by, 
for instance, bringing pairs of spheres together - mov- 
ing two spheres closer requires no additional energy, but 
by reducing the number of nearest neighbors the overall 
energy is reduced. To wit, the same physics drives the 
formation of a multiply-occupied crystal in the penetra- 
ble sphere model (cr = 0) [12(, while the ground state 
of a generalized exponential repulsion is a multiply- 
occupied fee lattice with a lattice constant independent 
of the average density. To explore the HCSS ensemble 
we introduce a lattice model with occupation rii = 0, 1 
at each site to enforce the hard core repulsion. The re- 
maining soft shoulder is characterized by Vss (r) and the 
density is set through the chemical potential /i in the 
Hamiltonian: 
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hard cores [llj. We develop a self-consistent field theory 



Note that the sum is not over nearest neighbors but over 
all sites - the range of the interaction is encoded in Vss- 
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FIG. 1: Theoretical phase diagram and simulation snapshots 
of some solid and liquid modulated phases in two dimensions. 
Shown are only the solid (M s ) and liquid (Mi) micelle phases, 
solid (L s ) and liquid (Li) lamellar phases, and the solid inverse 
micelle phase (IM S ). Here, a s ja = 5 where a s is the soft- 
shoulder radius and a the hard-core radius. The potential 
height is e, the pressure is p, and f3 = l/ksT. The dark circles 
are the particle hard cores and the diffuse circles represent the 
particle soft-shoulder. 

This model, with a square shoulder, was used to study 
electron liquids in weak magnetic fields . There it was 
a "toy" model of the true interactions, while here we 
study its consequences for generic V^ss |lig • To develop a 
mean-field theory, it behooves us to rewrite the partition 
function in terms of a continuous field <f> through the 
Hubbard- Stratonovich transformation [l6| 

z = j2jmcJ-^ J -0-j(k)v s - s \h)ci>(-h) 

+ X>j(#j+/W (3) 
j ) 

where we have mixed the continuum and discretum for 
notational ease and d is the dimension of space. Here 
Pa = ld&~ d is the density of the lattice where we have 
chosen the lattice constant to equal the hard sphere di- 
ameter <7 and 7d is a lattice and dimensionally dependent 
geometrical factor. Summing over rij = 0, 1 results in an 
interacting Fermi-Dirac-like free energy. The average site 
occupation is (n z ) = /^dlnS/d// = Ae l ^/(1 + Ae*) 
where A = e"^ is the fugacity. This is related to the 
number density profile through pi — pa{ni). Solving 



the mean-field equations, we find that V S g (k)0(k) = 
i(3p(k)/pQ and the mean-field estimate for the fixed p 
(not fixed p) free energy, obtained via Legendre trans- 
form, is the standard density functional form: 

PF= ^2 / (^P(k)fe(k)p(-k)+ J d^[p\n(p/p ) 

+ (p -p)\n(l-p/p )] . (4) 

Following Likos et al. we consider the stability 

of fluctuations around a uniform density. Working to 
quadratic order in fluctuations about a constant p we 
find an instability when both Vjg^k) < and po < 
p/po (1 — p/po) /3|Vss(k)|. At low volume fractions, when 
p/ Po ~ * anc i the hard cores rarely overlap, we recover 
the result in Ref . ^lj ■ The presence of the hard core not 
only alters the stability criterion, but introduces a new 
length scale a which competes with the overlap of the 
shoulders at a s . 

Though our lattice model provides a modicum of in- 
sight, it fails to distinguish ordered from disordered states 
in the clusters. Thus, we have analyzed the HCSS sys- 
tem off lattice in two dimensions by constructing a mean- 
field model which captures the salient features of the self- 
assembly. We consider the interaction of clusters of fixed 
shape in various periodic arrangements. The interaction 
between the clusters arises from the summed interactions 
of the spheres in each cluster, rotationally averaged so 
that we can ignore any spatial or orientational correla- 
tions between spheres in different clusters. It is then 
straightforward to calculate the average overlap energy 
for a given aggregate morphology. Moreover, for each 
morphology, we consider either fluid or crystalline order 
in the clusters; the entropy is estimated using an empir- 
ical equation of state for the hard-disk fluid or the 
cellular theory of the hard-disk crystal ^| . The result is 
a closed, albeit cumbersome, analytic form as a function 
of temperature, density, and the two structural parame- 
ters, cluster size and lattice spacing (the latter are fixed 
via Lagrange multipliers). An example of the phase di- 
agram of the aggregate phases, which are stable at low 
enough temperatures /3e > 1, is shown in Fig. ^ 

In units of crj, the wavenumber at which ordering 
occurs only depends weakly on the dimensionless pres- 
sure p&g/e and is, interestingly, roughly independent of 
the particular phase, be it micellar, lamellar, or inverted 
micellar, as shown in Fig. [21 This behavior recapitu- 
lates the results for a ge neralized exponential interaction 
without hard cores [13|. Taken together with the obser- 
vation that the qualitative features of the phase diagram 
are independent of (J s /cr, this confirms our notion that 
we can think of the clusters effectively as larger inter- 
acting particles with the internal structure only chang- 
ing the effective cluster-cluster interactions. The relative 
unimportance of the hard-core interaction arises when 
the clusters are large and an additional sphere can be 
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FIG. 2: Wavelength of principal mode vs. pressure for a s /a = 
5 and /3e = 2. The dashed line shows the prediction from the 
lattice, mean-field theory, L/a s ~ 1.22, the solid line shows 
the results of the mean-field, off-lattice calculation with the 
boundaries of the respective phases labeled below; and the 
filled circles show the results from the simulations with the 
phases labeled above. The phase morphologies are labeled as 
in Fig. Q The lattice constant and the wavelength differ by 
a multiplicative factor of y/3/2 in the hexagonal phases, as 
usual. 

accommodated without requiring a large rearrangement 
of the internal structure. Moreover, we find that in the 
pressure and temperature ranges we consider, the micel- 
lar and lamellar phases can have both fluid and crys- 
talline intra-aggregate order. In contrast, the only sta- 
ble inverted micellar phase is crystalline. Our analysis 
corroborates the results of Monte Carlo (MC) simula- 
tions of the HCSS system at fixed N, p, and T Q. 
At low temperatures, the MC simulations reveal a series 
of first-order transitions, with decreasing pressure, from 
crystalline inverse micelles, to crystalline lamellae, to liq- 
uid lamellae, to crystalline micelles, to liquid micelles, to 
isotropic liquid, in qualitative agreement with the phases 
and phase diagram of our mean-field treatment. Transi- 
tions between distinct structural variants of each type of 
phase (e.g., between different crystalline lamellar phases) 
are also observed. A more detailed discussion of the sim- 
ulation results appears elsewhere Q . 

To understand the unwavering magnitude of the 
wavevector, we first consider the lamellar phase in the 
limit where the boundaries between the occupied and 
unoccupied regions are sharp, in the spirit of the strong 
segregation limit of diblock copolymers || . Let £ be the 
width of the lamellae and a be width of the gaps between 
them so that the centers of the lamellae are L = £ + a 
apart. If £ < er s , the energy per sphere arising from 
other spheres in the same lamella is -Emtra ~ ep c g£a s 
where p c g = p(l + a)/£ is the density of spheres in the 
occupied region. Thus -E; n tra ~ ep(£ + a)o~s] a fixed lat- 
tice constant L amounts to fixing the energy per particle 
arising from intra-aggregate interactions. Qualitatively, 
the inter-lamellar repulsion decreases both with a and 



t. as a grows for fixed £, fewer particles in neighboring 
lamellae interact, and, as £ grows for fixed a, the fraction 
of the particles at the edge of each lamella which overlap 
more strongly with the neighboring lamella is reduced. 
To lowest order, our off-lattice mean-field analysis finds 
Winter ~ 0.75epy/a^(a s - a) b / 2 {£ + a)/£ 2 , and the mini- 
mum of -Einter + -Eintra (subject to the condition that p e g 
should not exceed the close-packed density 2/v / 3e>' 2 ) in- 
deed lies on the line £ + a ss l.2a s . We also find that the 
equilibrium value of £ is as small as allowed by this con- 
dition, as if the lamellae were being compactified by an 
effective surface tension arising from the inter-lamellar 
repulsion. 

Returning to the lattice model, we evaluate the free en- 
ergy of modulated phases; wc need only consider Fourier 
modes which belong to the reciprocal lattice of the pe- 
riodic structure (including the k = mode). Fixing the 
zero mode of the density, we find that those wavevec- 
tors k for which Vss(k) < are unstable, and that the 
strongest instability will be at the minimum (i.e., most 
negative) value of Vss- In the case of the soft shoulder 
in two dimensions, Vss(k) = eo~ s Ji(ko~ s )/k, which has a 
minimum at k*a s « 27r/1.22 in remarkably good agree- 
ment with the off-lattice mean-field results and the Monte 
Carlo simulations (Fig. [3) . 

When is it appropriate to assume the strong-segrega- 
tion regime and an infinitely narrow boundary between 
occupied and unoccupied sites? To assess the accuracy 
of this approximation, we return to the mean-field treat- 
ment of the lattice model at fixed chemical potential. 
In the low temperature regime, we expect the bound- 
aries to be infinitely narrow. In this limit, we can use a 
Sommerfeld-like expansion to compute p(r) for a stripe 
configuration to lowest order in A, the width of aggre- 
gate boundary. As is the case for the square shoulder, 
we assume that |Vss(k)| decays sufficiently rapidly with 
|k| and that the potential, </>(r), is well-approximated 
by 4>(y) ~ cj)Q + <fi*cos(k*x). These coefficients are re- 
lated by the mean-field equation to appropriate Fourier 
modes of the density. At sufficiently low temperature 
the density is a step function, and the relevant Fourier 
modes are p(k = 0)/V to t = Pok* £/2tt and p(k*)/V tot = 
po s'm(k*£/2)/ir, where Vtot is the total volume (note that 
(j)o and 4>* are not Fourier modes). Simultaneously, the 
nonlinear relation between p and (f> gives an additional 
constraint on the position of the interface, x = £/2: 
—ftp, = i<fio + i(p* cos(k*£/2). Combining this with the 
mean-field equations, we eliminate 4>o and (f>* to find a 
transcendental equation for £: 

V* k*£ 

— B mk*l = V —- t i, (5) 

where V = V ss (0)/ Po and V* = \V ss (k*)\/ Po . When 
we have a soft shoulder, Vq > V* and k*£/(2n) ~ p/V . 
Temperature does not strongly affect the width of the ag- 
gregates, but it does alter the sharpness of inner domain 
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boundary. The width of the interface is measured from 
the slope of the density at the interface, and is given by 
A -1 = |p _1 (dp/da;)| at x = i/2. To lowest order, we find 
fc*A/(27r) ~ k B T/[2V*sm 2 {Trp/p cS )}. Therefore, we can 
see that the aggregates "melt" from the boundaries as 
the temperature is raised due to the fact that the depth 
of the self-consistent potential is proportional to V* . We 
deduce a rough estimate of the stripe melting tempera- 
ture, T m , from the condition that k* A/2tt < 1 for ordered 
structures, or /csT m ~ V* sin 2 (7rp/ p e g). From this we see 
that the most stable structures occur at half-filling, when 
p c f[ = 2p, consistent with our instability criteria. These 
calculations can be repeated for the other morphologies 
by utilizing more modes with |k| = k*. 

Our lattice model can be further embellished by in- 
cluding a short-range attraction in addition to a long- 
range repulsion 0|. Though the details would be dif- 
ferent, it is straightforward to see that our conclusion of 
an ordering instability would remain and that our quali- 
tative arguments would still hold. Clearly, this aggrega- 
tion phenomenon is encoded in the Fourier transform of 
the soft part of the pair potential, and we expect that 
self-assembly within a broad class of repulsive systems 
can be described by this model. This analysis also sug- 
gests that the morphological similarities between repul- 
sive self-assembly and clustering due to com pe ting short- 
range attraction and long-range repulsions |2fj are more 
than superficial: they reflect the common structure of 
the Fourier transform of the interaction potential. 

In addition, we note that the lattice model, by virtue 
of its likeness to the Ising model, has a certain symmetry 
under toggling occupied and unoccupied sites (at half fill- 
ing). This gives us a way to understand the existence of 
the inverse phases, and suggests that this symmetry is ap- 
proximately preserved even in the off-lattice model. Fi- 
nally, we note that the observed phase sequence is analo- 
gous to the ordered phases of diblock copolymers [ilj and 
amphiphiles in water , despite the symmetry of the in- 
teraction. Near the critical point of the molten copolymer 
mean-field theory, quite generic Landau-theory consider- 
ations lead to the diversity of modulated structures in 
equilibrium j^. We expect the structure of the mean- 
field theory for the present colloidal model to be very 
similar to the copolymer theory, at least near the mean- 
field critical point. This suggests that, in three dimen- 
sions, HCSS systems can be expected to form spherical, 
cylindrical, and lamellar aggregates (and their inverses); 
indeed, such phases have been observed in MC simula- 
tions of the three-dimensional HCSS system [l4|. Diblock 
melts also exhibit a bicontinuous gyroid phase, and it 
would be intriguing to identify an analogue for repulsive, 
self-assembling particles. 
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